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ABSTRACT 

We analyze the linear, 3D response to tidal forcing of a disk that is thin and thermally stratified in 
the direction normal to the disk plane. We model the vertical disk structure locally as a polytrope which 
represents a disk of high optical depth. We solve the 3D gas-dynamic equations semi-analytically in the 
neighborhood of a Lindblad resonance. These solutions match asymptotically on to those valid away 
from resonances (previously obtained by Korycansky & Pringle 1995) and provide solutions valid at all 
radii r. We obtain the following results. 1) A variety of waves are launched at resonance, including r 
modes and g modes. However, the f mode carries more than 95% of the torque exerted at the resonance. 
2) These 3D waves collectively transport exactly the amount of angular momentum predicted by the 
standard 2D resonant torque formula. 3) Near resonance, the f mode behaves compressibly and occupies 
the full vertical extent of the disk. Away from resonance, the f mode behaves incompressibly, becomes 
confined near the surface of the disk, and, in the absence of other dissipation mechanisms, damps via 
shocks. In general, the radial length scale for this process is roughly r^/m (for resonant radius tl and 
azimuthal tidal forcing wavenumber to), independent of the disk thickness H. This wave channeling 
process is due to the variations of physical quantities in r and is not due to wave refraction. 4) However, 
the inwardly propagating f mode launched from an m = 2 inner Lindblad resonance experiences relatively 
minor channeling (accompanied by about a factor of 5 increase in nonlinearity) , all the way to the radial 
center of the disk. 

We conclude that for binary stars, tidally generated waves at Lindblad resonances in highly optically 
thick circumbinary disks are subject to strong nonlinear damping by the channeling mechanism, while 
those in circumstellar accretion disks are subject to weaker nonlinear effects. We also apply our results 
to waves excited by young planets for which to « r/H and conclude that the waves are damped on the 
scale of a few H. 

Subject headings: accretion: accretion disks — binaries: close — planetary systems — solar system: 
formation — stars: pre-main-sequence 
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1. INTRODUCTION 

Gaseous disks are found in many types of binary star 
systems, including cataclysmic variables (CVs) and pre- 
main-sequence stars, and young planetary systems. The 
orbiting objects (stars or planets) exert tidal forces on 
these disks, which generally act merely to distort the disks 
from an axisymmetric form. However, at special locations 
in disks where resonances occur, the tidal forces generate 
waves that transport energy and angular momentum. As a 
result, a resonant torque is exerted by the system objects. 
The orbital evolution of the perturbing objects sometimes 
depends on the strength of these torques (Goldreich & 
Tremaine 1980, hereafter GT80; Lin & Papaloizou 1993; 
Lubow & Artymowicz 1996). The waves transfer their an- 
gular momentum and energy to the disk in the regions of 
space where they damp, and this in turn affects the evo- 
lution of the disk. For example, gaps could be created in 
disks in regions where the waves damp. 

This paper concentrates on Lindblad resonances (LRs), 
which are due to horizontal forcing (along the disk plane). 
We assume throughout that the disk is coplanar with 
the orbit of the system objects. A 2D, linear theory for 



resonant tidal torques and associated wave propagation 
was developed by Goldreich & Tremaine (1979, hereafter 
GT79). This theory provides an explicit formula for the 
Lindblad resonant torque. The 2D theory considers the 
disk to have only radial and azimuthal extent and ignores 
effects over its vertical extent (perpendicular to the orbit 
plane) . Important progress has been made in the study of 
2D nonlinear waves in disks (Shu, Yuan, & Lissauer 1985; 
Yuan & Cassen 1994). The torque in the nonlinear case 
was found to be within a few percent of that predicted by 
the 2D linear formula. The nonlinearities produce highly 
spiked density profiles which increase the level of dissipa- 
tion present in a viscous disk. Radiative damping of linear 
waves can be important particularly when the disk is warm 
(Cassen & Woolum 1996). Turbulent viscosity in the disk 
provides another means of wave damping. 

The 2D treatment is valid if both the vertical structure 
of the disk and its thermodynamic response are locally 
isothermal. Under such circumstances, tidal forcing will 
generate a 2D wave in a 3D disk. The wave front remains 
perpendicular to the disk plane at all heights, as the wave 
propagates radially in the disk. However, this 2D wave is 
highly singular in that it does not exist in a disk with a ver- 
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tical temperature variation (Lin, Papaloizou, & Savonijc 
1990a, hereafter LPS; Lubow & Pringle 1993, hereafter 
LP). Furthermore, a vertically isothermal structure is not 
realistic for many important classes of gaseous disks, such 
as accretion disks in CVs, circumstellar and circumbinary 
disks of YSOs, and protoplanetary disks. Such disks of- 
ten have optical depths much greater than unity and can 
be expected to have substantial vertical temperature varia- 
tions, if they have an internal heat source such as turbulent 
dissipation. 

We demonstrate in this paper that 2D tidal forcing, 
caused by LRs, excites 3D waves in a thermally stratified 
disk. However, another class of resonances exist due to 3D 
effects. These resonances are due to vertical, tidal forcing 
by a coplanar perturber of a disk with nonzero thickness 
(Lubow 1981). The vertical resonances also generate hori- 
zontally propagating waves. Although intrinsically weaker 
than the LRs, the vertical resonances may be of impor- 
tance in close binary star systems. 

Some investigations of 3D effects have been carried out 
using numerical simulations (Lin, Papaloizou, & Savonije 
1990a, b). Such approaches can be used to explore a lim- 
ited range of physical parameter space. However, recent 
progress has been made in obtaining semi-analytic solu- 
tions for waves in 3D disks (LP; Korycansky & Pringle 
1995, hereafter KP). The aim of this paper is to extend 
that approach to understand 3D wave generation at LRs 
and the subsequent wave propagation. We determine the 
linear response of a thin (but nonzero thickness) disk. We 
model the disk locally as a polytrope in the vertical direc- 
tion, which is valid for a disk of very high optical depth, 
and we ignore the effects of atmospheric layers. We aim 
to understand which modes of a thermally stratified disk 
are excited at LRs and how much torque is carried by such 
waves. By studying the properties of linear wave propa- 
gation, we can also understand where nonlinearity sets in 
that will likely lead to shocks and subsequent wave dissi- 
pation. 

The outline of this paper is as follows. In §2, we review 
the properties of disk modes. In §3 and §4, we derive and 
solve the equations for waves generated at LRs. In §5, 
we compute the total torque carried by these waves and 
determine which modes are excited. In §6, we summarize 
properties of the f mode of the disk and present our nu- 
merical results. In §7, we discuss the application of our 
results to binary and protoplanetary systems. §8 contains 
a summary. 

2. REVIEW OF THREE-DIMENSIONAL WAVES IN A THIN 

DISK 

2.1. Free Waves 

In a thin disk, there is a separation of scales between 
the horizontal and vertical directions which implies that 
the waves take a WKB form in the radial direction, ex- 
cept in the neighborhood of resonances, which are also 
turning points for the waves. The equations governing ax- 
isymmetric WKB waves in a vertically isothermal disk are 
described in LP within the approximation of the shear- 
ing sheet, in which the radial derivatives of all base-state 
quantities other than the angular velocity are neglected. 
The case of a vertically polytropic disk was subsequently 
considered by KP. 



The problem reduces to a second-order system of or- 
dinary differential equations in the vertical coordinate z, 
which, together with appropriate boundary conditions, 
constitutes an eigenvalue problem for the frequency ui of 
the mode. The radial wavenumber k(r) appears as a pa- 
rameter, so that a local dispersion relation uj(k;r) is de- 
fined. Moreover, the same dispersion relation applies to 
non-axisymmetric waves if the frequency eigenvalue lo is 
replaced by the intrinsic frequency Cj = to — mfl, provided 
that the azimuthal wavenumber m > satisfies m <C r/H. 

The solutions of the eigenvalue problem can be classified 
by analogy with the theory of stellar oscillations (e.g. Cox 
1980). The name of each class of mode refers to the domi- 
nant restoring force that is involved in its dynamics. Thus 
p modes are acoustic modes, which are inherently com- 
pressible and propagate by pressure forces. The g modes 
are gravity modes, which rely on buoyancy forces resulting 
from an entropy gradient in the vertical direction. There 
are also r modes, which do not exist in a non-rotating star, 
and which propagate by inertial forces. The p, g, and r 
modes can be further classified according to their symme- 
try about the mid-plane and the number of nodes in the 
pressure perturbation (or radial velocity perturbation) in 
z > 0. Finally, the f mode is the fundamental mode of 
oscillation, and has unique characteristics described in de- 
tail in §6.1 below. In fact, a disk has two f modes, one of 
each symmetry, because it has two surfaces. 

Some of these modes are lost in special cases. In an in- 
compressible disk, there are no p modes, although all other 
modes survive. The g modes do not exist in a disk that 
is adiabatically stratified in the vertical direction. Simi- 
larly, there are no r modes in a disk that has zero epicyclic 
frequency. In an isothermal disk with no surface, the g 
modes are lost and the f mode of even symmetry takes a 
special form, known otherwise as the Lamb mode or two- 
dimensional mode. (There has been some confusion in ter- 
minology. LP used the naming convention that r modes 
were called g modes.) 

For a monochromatic wave or a normal mode of the 
disk, cj is independent of r and the dispersion relation can 
be followed continuously to determine the wavenumber at 
any radius. Typically, the wave occupies a finite interval 
in radius, which is bounded either by turning points, at 
which k goes to zero, or by the edges of the disk. More- 
over, the continuous variation of the amplitude of the wave 
with r can be determined by appealing to the conservation 
of either energy or angular momentum wave action. 

In many cases of interest, the dispersion relation is fol- 
lowed into a limit in which the dimensionlcss wavenumber 
kH becomes large. The behavior of the modes in this limit 
was described by Ogilvie (1998). Provided that the disk 
has a surface, the f, p, and g modes become confined in a 
layer of characteristic thickness k~ x near the surface, and 
their frequencies scale proportionally to fc 1//2 . In contrast, 
the r modes become confined in a layer of characteristic 
thickness fc -1 / 2 ^ 1 / 2 centered on the mid-plane, and their 
frequencies scale proportionally to fc~ 1//2 , unless the disk 
is marginally stable to convection. None of the modes oc- 
cupies the full vertical extent of the disk in this limit. 

It is useful to note the following properties of waves in 
a Keplerian disk. The f and g modes have \Co\ > and 
propagate both inside the inner LR and outside the outer 
LR. The r modes have \Cj\ < Q and propagate in the re- 
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gion between the two LRs, where the corotation resonance 
is also located. An exception is the 'tilt' mode, consisting 
of the union of the f and ri modes of odd symmetry, which 
propagates at all radii. The p modes have \lo\ > cj p , where 
w p > ft is a frequency which increases with the order of 
the mode. Their resonances are separated from the Lind- 
blad resonances by a distance which also increases with 
the order of the mode. These properties are summarized 
in Figure 1. 

2.2. Asymptotic Matching 

One of the goals of this paper is to determine to level 
of excitation of the free modes described above, due to 
tidal forcing at Lindblad resonances. To accomplish this, 
we determine in §4 the solutions for the tidally generated 
waves in a small region close to the Lindblad resonance (of 
order (i? 2 ?-) 1 / 3 in radius), where the wave excitation oc- 
curs. In this region of space, a locally valid approximation 
is used, in which the quantity k 2 — lu 2 , which has a zero 
at the resonance, is treated as a linear function of radius. 
Away from the resonance (over most of space) , a different 
approximation applies, which is a WKB approximation in 
the radial direction. Radial WKB free wave solutions of a 
vertically polytropic disk were obtained by KP. A solution 
that is uniformly valid for all r (asymptotically accurate 
in all space in the limit H <C r) is obtained by the method 
of matched asymptotic expansions (e.g. van Dyke 1975). 
In this context, the global wave solutions are determined 
by matching the outer limit (limit going away from reso- 
nance) of the driven Lindblad resonance solutions to the 
inner limit (limit approaching resonance) of the KP free 
waves. In practice, the level of excitation of each KP mode 
is easily determined from the amount of angular momen- 
tum generated for that mode in the region near resonance. 

3. MATHEMATICAL FORMULATION 

3.1. _Ba.se state 

Consider a cylindrical coordinate system (r, <f>, z) such 
that z — coincides with the disk mid-plane and the di- 
rection of increasing <p is the direction of disk rotation. 
The unperturbed disk has angular velocity £l(r), density 
p(r, z), and pressure p(r, z). We assume that the disk is 
thin and ignore effects of viscosity and self-gravity. As in 
KP, we adopt a polytropic disk structure in z locally in r 
so that 

p(r,z)= Po (r)(l-z 2 /H 2 y (1) 



and 



p(r, z)= Po (r)(l-z z /H<) 



2\s+l 



(2) 



where pa(r) and po(r) are arbitrary smooth functions of 
r and s is the polytropic index. The disk half-thickness 
H(r) satisfies 



paper and are generally ignored. The polytropic struc- 
ture provides a mathematically convenient description of 
a thermally stratified disk of high optical depth. However, 
the features of the results obtained in subsequent sections 
do not depend on the thermally stratified disk being poly- 
tropic. 

3.2. Dynamical equations 

Let the local Eulerian perturbations of disk velocity in 
cylindrical coordinates be represented by (it, v, w), that of 
pressure by p' , and that of density by p'. We represent 
physical wave quantities such as the pressure perturbation 
by p'(r, <j>, z, t) — p'(r, z) exp [im(<f> — f2 p f )], where positive 
integer m is the azimuthal wavenumber and 17 p is the an- 
gular pattern speed of the wave. Consider the disk to be 
subject to a tidal potential component < I > m (r, 0, z, t) of the 
above form. In the case that the system objects are in 
circular orbits, the pattern speed is equal to the angular 
orbital speed of the objects. If the system has eccentricity, 
then additional components of the perturbing potential 
appear containing the pattern speeds that differ from the 
fundamental orbital frequency of the system (see GT80). 
Using the 3D shearing sheet equations (see LP and KP), 
we obtain the following equations for the disk perturba- 
tions. 

on , dr P' d<Pm ca\ 
— iu)u — 2llv-\ = — , (4) 

p ar 



iojv + 2Bu = <I> m , 

r 



luiw + g — 
P 



d zP ' 



= 0, 



(5) 
(6) 



iui— + wd z In p + d r u + d z w = 0, (7) 
P 



and 



p p 

iuj 7 — 

p p 



• ri). In [ ) = 0, (8) 



where 7 is the usual adiabatic exponent. These equations 
are respectively the radial, azimuthal, and vertical force 
equations, the mass conservation equation, and the en- 
ergy equation for adiabatic perturbations. The frequency 
uj is defined as 



d)(r) = m [Q p — 0(t")] 



(9) 



Quantity B{r) is the usual Oort constant equal to ^ + §^r 
and g is the vertical gravity in the disk equal to g(r, z) = 
n 2 (r)z. The azimuthal pressure force imp 1 jr was ignored 
in equation (||) and the azimuthal velocity term imv jr was 
ignored in equation (^). This approximation is justified for 
a thin-disk limit for which 



H 2 (r) 



2(s + l)p (r) 
n 2 (r)p (r) ' 



(3) 



The thin-disk approximation used to derive the above re- 
quires H <C r. The base-state angular velocity is slightly 
modified from a Keplerian law by the radial pressure forces 
and it varies slightly in z within the disk. However, these 
effects are unimportant for the issues described in this 



H <C r/m. 



(10) 



The forcing term in the vertical force equation (g|) was 
ignored. This approximation can be justified by the thin- 
disk condition (fL(t) everywhere except at vertical resonance 
points (Lubow 1981). We concentrate on wave properties 
near Lindblad resonances (horizontal resonances) in this 
paper and thus this approximation is justified. 
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4. EQUATIONS NEAR LRS 
By combining equations (|J) and (|s|), we obtain 



u(lu 2 



drP' 



dG> 

dr 



2imfl 



-$ r . 



(11) 



where k is the epicyclic frequency such that k 2 = ABil. 
Notice that equations (§|)-(||), together with equation (fLl|), 
form a complete set, i.e., v has been eliminated. Once u 
has been determined, v can be obtained from equation (||). 

LRs occur at radii rx where Co 1 = k 2 . Away from such 
locations, the above equations are equivalent to the WKB 
equations of KP. Near an LR, the WKB approximation 
breaks down in equation ([ll]) , but the solutions can be rep- 
resented in terms of Airy functions in the neighborhood of 
the turning point. To treat this region, we follow standard 
techniques (e.g. GT79) and expand in the neighborhood 
of the resonance 

Vx, (12) 



2 -2 
K — U) 



where x = (r — ?"l)/Vl, and D = rd(K 2 — uj 2 )/dr is eval- 
uated at r — ix. We then approximate the radial force 
equation by 



uDx 



d r p' 



dr 



2im£l 



(13) 



4.1. Free waves 



The equations for free waves near LRs are obtained by 
setting $ m = and combining equations and ([l3|). 

We seek solutions of the following form. 



u(x, z) 
v(x, z) 
w(x, z) 
p'(x,z) 
p'{x,z) 



Ai{qx)u(z), 
iAi(qx)v(z), 
qAi' (qx)w(z), 
iqAi'(qx)p%z), 
iqAi'(qx)p'(z), 



(14) 



where Ai is the Airy function given by equation (10.4.1) 
of Abramowitz & Stegun (1965) and Ai' is the derivative 
with respect to its argument. The parameter q is a dimen- 
sionless constant that will be determined as an eigenvalue. 
(Solutions involving the complementary Airy function Bi 
have been discarded because, on one side of resonance, 
they approach infinity with distance from the resonance.) 

Using the property that Ai"{x) = xAi(x), we obtain the 
following free-wave equations for u, v, w,p' and p' . 



" 3 

pVu + — p' = 0, 

r 



IB. 



-u, 



ujpw - gp 



dp' 



= 0, 



„ -/ .dp 
ujp + w— + p 

dz \r 



dz 
u dw 
dz 



and 



P' P' 
w 7— 

\P P 



(15) 

(16) 
(17) 
(18) 

(19) 



where all coefficients are evaluated at r = r L . 



As can be seen from the above, the equations have been 
reduced to ID equations in z. We introduce dimensionlcss 
variables 



Z 

X(Z) 
F 

Q 
s 

W{Z) 



z/H, 

u;p'(z)/(n 3 H 2 p), 
{H/rf/Zq, 

v/n 2 , 

w(z)/(CLH). 



(20) 



Combining the above equations (^5|)-(p0|), we obtain 



dX 
~dZ 



and 



dW 
~dZ 



2aZ 
1-Z 2 



2(3Z 
1-Z 2 



X — \ F — 



2aZ 2 
1-Z 2 



W = 



w 



2/3 



1 - Z 2 



91 
s 



X = 0, 



(21) 



(22) 



where the constants a = s — (1 + s)/-f and j3 = (1 + s)/-f 
are proportional to the square of the disk's buoyancy fre- 
quency and the square of the inverse sound speed at mid- 
plane, respectively. We note that these equations can be 
derived from equations (13) and (14) of KP in an expan- 
sion about the point of the dispersion relation where k = 
and F 2 — k 2 . 



Equations (21) and ( |22| ) are integrated numerically with 
appropriate boundary conditions at Z — and Z — 1, as 
in KP. Both even and odd solutions about Z = can be 
obtained, but for the coplanar forcing problem at hand, 
only even solutions are relevant. Such solutions must then 
satisfy W-^O) = 0. Infinitely many solutions exist which 
differ in n, the number of vertical nodes of X (the pressure 
perturbation) for Z > 0. At Z = 1, the boundary condi- 
tion is that the Lagrangian pressure perturbation must 
vanish. Solutions are obtained by numerically integrating 
out of the regular singular point Z = 1 down to Z = 0. 
Only certain values of Q permit the mid-plane boundary 
condition W = to be satisfied. Once the value of Q 
is determined, the spatial structure of all physical quan- 
tities is known. Solutions with various spatial structures 
are obtained that behave as g modes, r modes, and the 
f mode. Each solution is characterized the type of mode 
(g, r, etc.) and by the value of n. There are a countably 
infinite number of solutions, which we label by index j. 

We note that the scalings for a thin disk are such that 
q is of order (r/iJ) 2 / 3 . The characteristic width of the 
resonance is therefore (7f 2 r) 1 / 3 , intermediate between H 
and r. In addition, the vertical component of the velocity 
is smaller than the horizontal components by a factor of 
order (H/r) 1 / 3 . 

4.2. Driven waves 

We return to the problem of wave generation due to 
tidal forcing. We consider solutions to the inhomogeneous 
equations near resonance, equations (||)-(||) an d (13). We 
consider the set of all homogeneous solutions obtained in 
the last section, each of which is characterized by the set 
of physical quantities qj, Uj(z), Vj(z), uij(z), p'Az), and 
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p'j(z). We construct possible solutions to the inhomoge- 
neous equations which are expanded in terms of the ho- 
mogeneous solutions in the following way. 



u(x, z) 
w(x, z) 
p'(x,z) 
p'[x,z) 



Ej a jUj(z) [Ai( qj x) + isjGi(qjx)} , 
J^ajqjWjiz) [Ai'{ qj x) + is 3 Gi' (qjx)] , 
i L,j a jQjPj( z ) [ Ai '(lj x ) + isjGi'iqjx)} , 
* Ej "j'hi'y-'- 1 [ Ai '(ljx) + isjGi'iqjx)} , 

(23) 

where Gi is the inhomogeneous Airy function defined in 
equation (10.4.55) of Abramowitz & Stegun (1965), Sj is 
either +1 or —1, and aj are the expansion coefficients. The 
linear combination of Airy functions, Ai(qjx) + isjGi(qjx) 
and Ai'(qjx) + iSjGi' , in the above is chosen to provide 
traveling waves in the limit qjX — > -co. The sign of Sj in 
the above determines the radial propagation direction of 
the wave (toward or away from the resonance) which we 
leave as an open issue for now for each j. 

We use the fact that Gi satisfies Gi"(x) — xGi(x) = 
—1/tt, and substitute the above forms into inhomogeneous 
equations d)-(i) and @. We find that equations (§)-(§) 
are automatically satisfied and do not in any way constrain 
the expansion coefficients. The constraint comes from the 
force equation (jl3|), which is satisfied provided that 



3 



2 P'j(z) 
p(z)- 



where forcing term is defined by 



d<f>, 



2m,n 



dr 



-<E> 



(24) 



(25) 



evaluated at r = r^. In equation (p4|), the left-hand side 
is independent of z. In the case of an isothermal disk 
with adiabatic index equal to unity (7 ~ 1), the ratio 
p'j/p is independent of z (and proportional to the sound 
speed squared) for the 2D mode, which is then the proper 
solution in that case only. In general, for disks with non- 
isothermal vertical structures or vertically isothermal disks 
having 7^1, the quantity p'J p varies with z for all modes 
in the disk. So some suitable combination of solutions is 
required. 

To solve equation (|24|), we seek an appropriate inner 
product that allows us to invert that equation to solve for 



aj. This can be derived by noting that equations (|15|)-([19|) 
can be combined to give a single equation for u in the form 



d / du 
dz \ dz 



+ fiu - 



Xpu 

H 2 n 2 



0. 



where 



N 2 



q *H 2 n 2 

r 2 V 



(26) 



(27) 



(28) 



is a dimensionless eigenvalue, and N is the vertical buoy- 
ancy frequency in the disk. (The detailed form of ji is not 
illuminating.) Equation (|6() does not satisfy the condi- 
tions of the Sturm-Liouville theorem, because /1 is not of 
definite sign and both /1 and fi vanish at z — ±H . (We 
note in passing that the singular point where u) 2 = N 2 is 



only an apparent singularity, since both linearly indepen- 
dent solutions are in fact regular there.) However, it can 
still be demonstrated that the eigenvalues are real 

and the eigenfunctions {uj} orthogonal in the sense that 



pu*Uk dz = 



(29) 



for any two modes with distinct eigenvalues. (Here and 
below, integrals are taken over the full vertical extent of 
the disk.) This is a limiting case of a more general orthog- 
onality relation, 



p (u*u k + w*w k ) dz = 0, 



(30) 



which applies to the eigenfunctions away from resonances 
studied by KP. Moreover, it can be asserted that the eigen- 
functions form a complete set on the space of continuous 
functions on (— H, H), although, technically, the pointwise 
convergence of an eigenfunction expansion is non-uniform 
with respect to z at the end-points where the weight func- 
tion p vanishes. 

By rewriting equation ( |2~i| ) in the form 



^^s^-Ujiz), 



<lj 



(31) 



and using the orthogonality relation, we obtain the coeffi- 
cients as 



aj = — / pu s dz 



p\uj\ 2 dz. 



(32) 



5. DERIVATION OF THE TORQUE FORMULA 

We now proceed to show that the torque exerted by 
a companion object at a Lindblad resonance is precisely 
equal to the value derived under a two-dimensional approx- 
imation. The radial velocity perturbation corresponding 
to the collection of waves excited at the resonance is known 
from §4.2 to have an inner expansion 



u(x, z) — ajUj(z) [Ai(qjx) + isjGi(qjx)} 



(33) 



in the vicinity of the resonance, where Sj = ±1 is still to be 
determined, and a factor exp(— iuot + im(j)) is understood. 
From the asymptotic forms (Abramowitz & Stegun 1965) 



-00 : Ai(x) ± iGi(x) 
-00: Ai(x)±iGi(x) 
x exp -JVz 



±iTT- 1 x~ 1 , (34) 
±iir- 1 / 2 (~xy 1 / 4 
;(-x) 3 / 2 + ^]}, (35) 



it is clear that the component of the solution labeled j 
matches on to a WKB traveling wave in the region q^x < 
and on to a decaying disturbance in the region qjx > 0. 

For single WKB waves of the disk it can be shown that 
the flux of angular momentum wave action, averaged over 
i, integrated over <f> and z, and measured in the direction 
of increasing r, is 



p\u\ 2 dz. 



(36) 
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The density of angular momentum wave action, similarly 
averaged and integrated, is 



out the possibility that sgn(2?) = sgn(J)), and so a more 
general formula is 



dz. 



(37) 



T 



-sgn(w) 



W\ 



(44) 



and a straightforward application of first-order perturba- 
tion theory confirms that the ratio /A^ is equal to 
the radial group velocity 



1)k 



(u> 2 



p\u\ 2 dz / J p{\u\ 2 + \w\ 2 ) dz 

(38) 

which is the slope of the WKB dispersion relation. We 
therefore identify the direction of propagation of a WKB 
wave according to the sign of its radial group velocity. 

Consider first the outer limit x — > — oo of equation (33) 
for a component with qj > 0. The wavenumber, computed 
as the derivative of the phase, is 



k 



It follows that 



sgn(u g ) = Sj-sgn(P)sgn(w), 



(39) 



(40) 



and, since sgn(2?) = — sgn(a)) for a Keplerian disk (with 
the convention that m > 0), the outgoing wave that is 
required by causality is the one with Sj = +1. From equa- 
tions (IH), and (j36|), the angular momentum flux 
associated with the wave is 



F 



w ,^fd! /(i|ili| ..„. 



(41) 



although we emphasize that this is measured in the direc- 
tion of increasing r, while the wave propagates into r < r^. 
By a similar argument it can be shown that, for a com- 
ponent with qj < 0, the required wave has Sj = —1 and 
carries an angular momentum flux 



F 



(a) 



<- 2 mV\a^ 
Lb 2 q 2 



p\uj\ dz 



(42) 



into r > j-l. From equation (p2|), the torque exerted on 
the disk is therefore 



T 



r 2 mT> ^ |a,-| , ., , 

L4r / p\ u j\ dz 



■K 2 m^ 2 



V 



E 



puj dz 



p\v,j\ dz 



V 



(43) 



where El = J pdz is the surface density at the resonance, 
and we have made use of the completeness of the eigen- 
functions {iij}. This formula is in precise agreement with 
equation (A10) of GT79. This result is not surprising since 
the 2D torque formula holds under a variety of physical 
conditions (Goldreich & Tremaine 1982; Papaloizou & Lin 
1984). For non-Keplerian rotation curves, one cannot rule 



We now introduce a specific disk model which is used 
for all the numerical calculations in this paper. The disk is 
Keplerian and its vertical structure is polytropic, with in- 
dex s = 3. The adiabatic exponent is taken to be 7 = 5/3. 
For this model we have computed the fraction 



fj 



puj dz 



El / p\v,j\ 2 dz 



(45) 



of the torque carried by each mode. The results for the f 
mode and the first three g and r modes are presented in 
Table 1, together with the corresponding eigenvalues Xj 
(defined by eq. (^8|)). It is clear that almost all the angu- 
lar momentum is transported by the f mode, although the 
r and g modes are also weakly excited. These first seven 
modes together account for 99.93% of the torque. For ob- 
vious reasons we focus on the propagation of the f mode 
for the remainder of this paper. However, it is worth not- 
ing that the r modes propagate away from the resonance 
in the opposite direction to the f and g modes. The results 
in Table 1 are independent of m. 

6. PROPAGATION OF MODES IN A KEPLERIAN DISK 
6.1. Basic properties of the f mode 

In this section, we summarize the properties of the sin- 
gle most important mode in this problem, the f e mode or 
fundamental mode of even symmetry. The eigenfunction 
of this mode, whether at the resonance or away from it, 
is uniquely characterized as follows. The pressure pertur- 
bation (or radial velocity perturbation) has no nodes in 
—H < z < H . The vertical velocity perturbation has only 
a single node at z = 0, as required by symmetry. This 
form of the eigenfunction allows a large overlap integral 
with the forcing potential which is independent of z, and 
explains why this mode is excited more than any other 
(see eq. @). _ _ 

The dispersion relation for this mode in our model 
disk (which is Keplerian and polytropic, with s = 3 and 
7 = 5/3) is shown in Figure 2a. Also shown is the asymp- 
totic form 

uj 2 - n 2 kH (46) 

derived by Ogilvie (1998), which is valid in the limit kH — > 
00. Since the surface gravity of the disk is g s = Q 2 H, the 
above equation is consistent with the well-known disper- 
sion relation for surface waves in a deep layer of water. 
A more accurate approximation, which takes into account 
the effects of rotation and stratification, is 



n 2 kH 



2 « 2 



(47) 



and this is shown in Figure 2b to provide an excellent ap- 
proximation even for values of kH as small as 3. 

The form of the eigenfunction in this limit also corre- 
sponds to that expected of a surface gravity mode. For def- 
initeness, we describe a wave traveling radially outwards 
in the disk, but inward traveling waves and standing waves 
are very similar in form. The wave is trapped in a layer of 
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Table 1 

Modes excited at a Lindblad resonance. 



Mode 




fi 


f 


6.16 


0.9678 


gi 


112.26 


0.0010 


g2 


236.47 


0.0003 


g3 


401.73 


0.0002 


I'l 


-18.98 


0.0277 


T2 


-98.32 


0.0018 


T3 


-237.42 


0.0005 



characteristic thickness k 1 near the surface of the disk, 
and in this layer the asymptotic forms 

u ~ A cos(kr — Lot + m(j>) exp(— kS), (48) 
w ~ A sin(fcr — tot + mcf)) exp(— kS) (49) 

are valid, where S = H — \z\. This velocity field is both 
irrotational and solenoidal (i.e. non-compressive) . The az- 
imuthal velocity is smaller by a factor O (ikH)^ 1 / 2 ) . In 
this limit, the even and odd f modes are indistinguishable 
since there is a loss of contact between the two surfaces of 
the disk (KP). 

For small values of kH, the mode departs from the be- 
havior described above. An approximation valid in the 
limit kH — > is 

to 2 ~ k 2 + k 2 H 2 n 2 /\, (50) 

where A > is the dimensionless f-mode eigenvalue de- 
fined by equation (|28|) . This form is reminiscent of the 
dispersion relation or = k 2 + c 2 k 2 for the two-dimensional 
mode in an isothermal disk (LP). While it is true that A 
depends on 7, implying that the f e mode is partly com- 
pressive in this limit, nevertheless A has a non-trivial limit 
when 7 — > 00. Therefore the mode behaves partly com- 
pressibly near the resonance, but only when the disk is 
compressible. The mode does not disappear in the incom- 
pressible limit. 

We emphasize how this description differs from ear- 
lier, less accurate theories of waves in disks. If the disk 
is described as a two-dimensional fluid, then all radially 
propagating modes are inherently compressive by virtue 
of the constraints on the fluid, and would normally be de- 
scribed as sound waves modified by rotation. In a three- 
dimensional disk, a two-dimensional mode of this form ex- 
ists only if the disk is vertically isothermal. In general, 
however, the principal mode excited at an LR in a three- 
dimensional disk is the f e mode, which is not inherently 
compressive. Moreover, this mode cannot be described 
using a vertical WKB approximation (cf. Vishniac & Di- 
amond 1989; Goodman 1993) because it has no pressure 
nodes in the vertical direction, being, if anything, evanes- 
cent rather than wavelike. 



6.2. Propagation and damping of the f mode 

6.2.1. Remarks on the nonlinear development of the f 

mode 

If the waves excited at a resonance were able to propa- 
gate radially without dissipation, they would reflect from 
the radial boundaries and establish a standing wave pat- 
tern. The net transfer of angular momentum and energy 
between the disk and the companion would then be zero. 
In order for the waves to have a significant dynamical ef- 
fect on the evolution of the system, they must deposit 
their angular momentum (either positive or negative) in 
the disk as a result of being damped. This can happen 
within linear theory if there is sufficient viscosity, but oth- 
erwise nonlinear dissipation is required. As noted by KP, 
the concentration of modes near the surface of a polytropic 
disk is likely to result in enhanced dissipation through the 
formation of shocks. 

It might be argued that, since the f mode is approxi- 
mately non-compressive, it would form shocks less easily 
and would therefore be less susceptible to this type of dis- 
sipation. However, little is known about the nonlinear de- 
velopment of surface gravity modes in a compressible fluid. 
In the case of an incompressible fluid, it is well known that 
when the amplitude of a surface gravity wave is increased, 
the distorted profile of the surface changes from a sinu- 
soidal shape to one with sharper crests, forming an angle 
approaching 120°. This happens when the velocity per- 
turbation is comparable to the phase velocity {gs/k) 1 / 2 
(e.g. Stokes 1880; Lamb 1932), and causes the crests to 
break. Although we consider it more probable that the 
nonlinear development in a compressible fluid will lead to 
shocks rather than breaking crests, it is still appropriate 
to normalize the solution described by equations ( |4"s| ) and 
(|49| ) according to 

/ Cl 2 H\ 

u ~ J\f ( — - — J cos(fcr — tot + mcf)) exp(— kS), (51) 

where M is a measure of the nonlinearity of the wave. In- 
deed, when this is substituted into the nonlinear equations, 
it can be seen that the nonlinear terms become compara- 
ble to the linear ones when |AA| exp(— kS) « 1. Therefore 
the wave can be considered to be marginally nonlinear at 
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the surface when \Af \ ~ 1, but, when |7V| ~ e, say, a sub- 
stantial fraction of the wave action can be considered to 
be subject to nonlinearity. 

6.2.2. Approximate, analytical treatment of radial 
propagation 

To be specific, we now consider the case of objects in 
a circular orbit. (The results for non-circular orbits that 
provide eccentric LRs do not differ in a substantial way 
from the analysis in this and the following subsection.) 
There is a Keplerian disk in which the angular velocity is 



1/2 



<>, [ - 

. r c . 



-3/2 



(52) 



where r c and £l c are the radius and angular velocity, re- 
spectively, of the corotation radius. The dimensionless 
intrinsic frequency of a mode with angular pattern speed 
fl c is 

3/2 

1 , (53) 

r r 



9. 



and the radii of the LRs are 



n, = 



m T 1 



2/3 



(54) 



where the upper and lower signs refer to the ILR and OLR, 
respectively. By combining equation (|53| ) with the dimen- 
sionless dispersion relation (Figure 2), we can derive the 
variation of the dimensionless wavenumber kH with r/r c . 
The results for various values of m are shown in Figure 3. 
It is clear that, after the wave has propagated over a radial 
distance comparable to r^/m from the resonance, kH has 
become sufficiently large that the wave can be described 
accurately as a surface gravity mode. This makes it pos- 
sible to determine the amplitude of the wave analytically. 

Let / denote the fraction of the torque carried by the 
f c mode. Then, from equation (]36|), the flux of angular 
momentum wave action associated with the wave is 



p\u\ 2 dz = ^ffT. 



(55) 



In the limit in which the wave can be described as a surface 
gravity mode, kH can be approximated as u) 2 /f2 2 , which 
itself can be considered large compared to unity, while the 
above integral becomes 



P\ 



u\ 2 dz 



2p Af 2 



n 2 H 



= r{s + i)M 2 (kHy 



26 
~H 

p n 2 H 



-2fc<5 



d6 



(56) 



If we assume that H oc r and po oc r p , the approximate 
radial variation of Af is given by 



AT 2 oc 



3/2 

-1 -1 



-i 2(s+3) 



p-2 



(57) 



As the wave approaches the radial center of the disk, 
TV 2 oc r p ~ 2 . It is quite plausible that p > 2 (implying that 



the surface density increases inwards faster than r^ 1 ), in 
which case the wave becomes less nonlinear as it travels in- 
wards. Otherwise, if the surface density increases inwards 
less rapidly, AT 2 increases inwards no faster than r -1 . For 
many cases of interest, there is therefore nothing singular 
about the wave approaching the radial center of the disk. 

As the wave propagates outwards into r S> r c , Af 2 oc 
r P+3s+7_ Thjg means that, for any reasonable density pro- 
file and polytropic index, the wave is very likely to become 
nonlinear as it propagates outwards. 

6.3. Numerical results 

We have computed numerically the radial propagation 
of the f c mode in a Keplerian disk using the method out- 
lined in the analytical treatment above, except that we do 
not apply the limiting analytic forms for surface gravity 
waves described after equation (55). The model disk has 
the standard parameters s = 3, 7 = 5/3, and also H oc r 



and po oc 



The results are presented as color images 



in Plates 1 and 2, for the cases m — 2 and m — 10, re- 
spectively. In each case, panels (a) and (b) refer to the 
inward propagation from the ILR, with radii < r < r^ 
shown. Panels (c) and (d) refer to the outward propaga- 
tion from the OLR, with radii r^ < r < 2r^ shown. In 
order to show the detail clearly, the vertical scale is greatly 
exaggerated. The two quantities plotted are the density of 
angular momentum wave action, averaged over time, 



and the RMS velocity perturbation, 



"RMS 



w 



2 )] 



1/2 



(58) 



(59) 



In panels (a) and (c), we plot the logarithm (to the base 



10) of \A^\, in units of \T\Q- L r^ l H^. In panels (b) 
and (d), we plot the logarithm of «RMs/(^-ff), m units of 
^11/2^-1/2^-1/2^-3/2^-1 The former quantity showg 

how the activity of the wave is distributed, in both the ver- 
tical and the radial directions. The latter quantity, apart 
from a normalization factor, is intended as a measure of the 
velocity amplitude of the wave relative to the local sound 
speed at the mid-plane. It is a more conservative measure 
of nonlinearity than Af 2 , and may be more appropriate in 
regions close to the resonance, where the f c mode behaves 
partly compressibly and may produce shocks more easily. 

If the disk extends continuously over the resonance, the 
torque can be calculated in terms of the standard formula, 
equation (43). In many cases of interest, however, the 
LR truncates the disk. In such cases, the disk edge par- 
tially overlaps with the resonance and the torque is often 
physically determined by the viscosity. The density at the 
resonance then adjusts so that its average value over the 
resonance width provides the needed torque in accordance 
with equation (J43|) . The meaning of El in the normal- 
ization of the velocity is the surface density expected on 
the basis of the smooth power-law density profile in the 
main body of the disk, rather than its actual value near 
the LR. In the case of a disk with an edge, the results in 
Plates 1 and 2 apply where the power-law density distri- 
bution holds in the disk. 
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Some care is required in interpreting Plates 1 and 2. 
They are plotted as if H/r were equal to 1/2 (for the ILR) 
or 1/4 (for the OLR), but this is done to allow a conve- 
nient visualization and has no intrinsic significance. The 
value of H/r is immaterial provided only that it is small; 
indeed the results become more accurate as H/r — > 0. 
The channeling of the wave occurs on a radial scale that 
is independent of H. There is a region of characteristic 
width (iJ 2 r) 1 / 3 around the resonance where, properly, the 
Airy functions should be used, but no attempt has been 
made to include this on the plots because it would require 
a definite choice to be made for the value of H/r. In fact, 
the singularity of the WKB solutions at the resonance is 
barely evident on the plots. 

For the case m = 2 (Plate 1), the inward propagation 
leads to a mild concentration of the wave action near the 
surfaces of the disk. However, the amplitude of the wave 
near the surface does not increase by large amounts away 
from the resonance and in fact tends to a constant value 
as the radial center of the disk is approached, in accor- 
dance with equation (|57|). The (dimensionless) amplitude 
increases by about a factor of 5 from the neighborhood of 
the LR to the disk center. With a larger value of p (more 
rapid density increase toward the center), the overall am- 
plitude increase would be less and the amplitude would 
decrease near the center. By contrast, a much more signif- 
icant channeling effect occurs in the outward propagation, 
and the amplitude near the surface increases strongly away 
from the resonance, also in accordance with equation ( p7| ) . 

For the case m = 10 (Plate 2), the channeling effect is 
very pronounced for both inward and outward propaga- 
tion. In the regions close to the resonances, there is an 
approximate symmetry between the two cases. However, 
it is again seen that the amplitude tends to a (large) con- 
stant value as r — > in the inward propagation, while the 
amplitude increases strongly as r — > oo in the outward 
propagation. 

7. APPLICATION TO BINARY AND PROTOPLANETARY 
SYSTEMS 

In this section, we explore the consequences of this work 
for some astrophysical situations. In any binary star sys- 
tem, up to three disks can occur. There may be two cir- 
cumstellar (CS) disks, each of which orbits about one the 
stars. In addition, there may be a circumbinary (CB) disk 
that orbits around the entire system. Powerful tidal forces 
present in binary systems clear large gaps in the disks, 
of order a few times the binary semi-major axis (Lin & 
Papaloizou 1993; Artymowicz & Lubow 1994, hereafter 
AL), which separate the disks from each other. As a con- 
sequence of the large size of the gaps, the disks experi- 
ence tidal forcing primarily from relatively low azimuthal 
wavenumbers m ~ 2. 

As discussed in the Introduction, there are several dissi- 
pation mechanisms that can act to damp the waves. This 
section contains some discussion of turbulent viscous dissi- 
pation of waves. The nature of the disk turbulence is quite 
uncertain. Additionally, the interaction of the f mode with 
turbulence is complicated by the fact that the largest tur- 
bulent eddies carry most of the viscosity, but they may 
not interact efficiently to damp the waves. Some simple- 
minded arguments suggest that the spatial damping rate 
in the radial direction is kj ~ ak(kH) b , where a is the 



usual parameterization of viscosity and 6 is a number of 
small absolute value that depends on the details of the 
turbulence. Such complications add to the uncertainty in 
k\, but we adopt kj ~ ak as a plausible estimate. 

In general, we expect wave channeling to amplify the 
wave substantially over a radial length scale r^/m for 
which kH ;> 1. A rough relative measure of the impor- 
tance of turbulent wave damping to wave channeling is 
kiTh/m ~ art,/(mH). For protostellar or protoplanetary 
disks, this ratio is less than unity, while for CV disks it is 
greater than unity, as is discussed below. 

One important example of CB disks occurs in young 
binary stars, where the remnant material involved in the 
star-formation process orbits around the binary. Typically 
these binaries have eccentricities e ~ 0.3 and the LRs are 
associated with the eccentric motions of the system (AL). 
An LR torque is responsible for truncating the CB disk 
at its inner edge. Such CB disks have been inferred to 
exist from the spectral energy distribution of some young 
binaries (Mathieu 1994). Based on our results in Plate 1, 
we expect that significant wave channeling would occur 
in a thermally stratified CB disk within a distance ~ tl/2 
from the resonance. We expect that strong nonlinear wave 
damping would occur in a CB disk, but exactly where 
this would occur depends on the level of nonlinearity of 
the wave near the LR. The 2D simulations of CB disks 
using smoothed-particle hydrodynamics (SPH) show that 
the wave is somewhat nonlinear near resonance, as is ap- 
parent from the prominence of the spiral arms near the 
disk's inner edge (AL). This effect provides some justifi- 
cation for the assumption that the wave damping occurs 
relatively close to the resonance (as occurs in 2D SPH 
simulations of CB disks, but possibly due to artificial vis- 
cosity). 

In the case of CS disks in young (eccentric) binary sys- 
tems, the eccentric Lindblad resonances can be important 
(GT80; AL). As seen in Plate 1, the wave channeling ef- 
fect is relatively weak, but of possible importance as the 
level of wave nonlinearity increases by about a factor of 5 
to the disk center. For low amplitude waves in the outer 
parts of the disk, the wave may well survive to the center 
of the disk, provided that other damping effects can be ne- 
glected. If the wave reflects off the central star, then the 
process of wave generation can become reversible, as the 
reflected wave returns to the LR. The net effect is a zero 
LR torque. Other wave damping effects, such as viscous 
wave damping or damping within the star, may play a role 
in preventing reversibility. 

In protostellar disks, the value of a has been estimated 
to be of order 10 -2 (e.g. Hartmann et al. 1998). The ex- 
tent of the turbulent decay of a wave over radius compara- 
ble to tl can be estimated by taking kH ~ m 2 , H/r ~ 0.1, 
and m ~ 2 for binary star forcing. The wave amplitude 
decays by roughly a factor exp (— kiri/) ~ 0.7. This plausi- 
bly provides some damping in a protostellar CS disk. For a 
CB disk, nonlinear damping by wave channeling will likely 
dominate, if the wave is slightly nonlinear near resonance. 

In addition, several classes of mass-transfer binaries (e.g. 
X-ray binaries and CVs) possess CS disks. These close bi- 
naries have nearly zero eccentricity, as a consequence of 
the tidal interactions between the system objects. Under 
these conditions, there are no eccentric LRs and the disk 
is sufficiently truncated that no (non-eccentric) LRs can 
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lie within the the disk. However, the closest LR, corre- 
sponding to to = 2, is sufficiently powerful to generate 
waves off-resonance as a result of the non-zero width of 
the resonance (see, e.g., Savonije, Papaloizou, & Lin 1994). 
Again, our results suggest that wave channeling is of pos- 
sible importance, but is not very powerful. In practice, 
the vertical resonances may play a more important role 
in such systems, since these resonances can lie within the 
disk (Lubow 1981). 

CV disks surround white dwarf stars and the value of 
a has been estimated to be of order 0.1 during outbursts, 
but could be much lower during quiescence (e.g. Can- 
nizzo 1993). The extent of turbulent radial decay rate 
over radius comparable to tl can be estimated by taking 
kH ~ to 2 , H/r ~ 0.02, and to ~ 2 for binary star forc- 
ing. The wave amplitude radial decay over a radius is of 
order 10 s . These estimates are crude, but suggest that 
the turbulent wave attenuation is severe (at least during 
outbursts), more severe than wave channeling, and would 
damp the f mode. 

In protoplanetary disks, the planets carve out relatively 
small gaps or possibly no gap at all. As a result, several 
high-TO LRs can lie within the disk. Under such condi- 
tions, it has been shown that the dominant torque comes 
from the LR having to w r/H (GT80; Ward 1986, 1988; 
Artymowicz 1993a, b). For typical protoplanetary disks, 
we then expect to <~ 10—20 to be most important. Our cal- 
culations have assumed that to -C r/H, so that we cannot 
apply our analysis throughout the resonant region. In ad- 
dition, these LRs have a vertical forcing contribution that 
is comparable to the horizontal forcing, and both of these 
vary somewhat with z (Ward 1988; Artymowicz 1993a). 
The present analysis assumed vertically constant horizon- 
tal forcing and ignored some azimuthal terms in the dy- 
namical equations, as discussed in §3.2. In principle, our 
methods could be extended to handle these complications. 

However, even in this case, some important conclusions 
can be drawn from the analysis in this paper. The region 
of resonant forcing in this case can be shown to be of order 
H in both the radial and vertical directions (rather than 
(i7 2 r) 1//3 for the radial and H for the vertical directions in 
the case to <C r/H). Furthermore, it can be shown that 
the waves for which to w r/H do approximately satisfy the 
equations of this paper at distances of order a few times 
H from the LR. This is because the radial wavenumbcr in- 
creases so rapidly with increasing distance from resonance 
that it quickly dominates the azimuthal wavenumber and 
makes our approximations valid. The f mode suffers wave 
channeling at radial distances from the resonance of order 
tl/to ~ H in this case. 

Therefore, we are led to the simple picture that the f 
mode is launched in protoplanetary disks from a region of 
size of order H. It undergoes wave channeling on a sim- 
ilar scale of a few H. As can be seen from Plate 2 for 
to = 10 (which we argue is approximately valid at dis- 
tances greater than about tl/5 from n, for H/r ~ 0.1), 
strong nonlinear damping can then occur with increasing 
distance from the LR. 

We expect that a full 3D analysis of the resonant re- 
gion would reveal that the f mode is still strongly excited 
in protoplanetary disks, compared with other modes. The 
reason is that the vertical component of the tidal forcing is 



not dominant and has a single zero at mid-plane. Further- 
more, the vertical variation of the horizontal forcing within 
the disk is not large, especially when the density weighting 
is taken into account (Ward 1986; Artymowicz 1993a). So 
the tidal forcing should preferentially drive modes with no 
or few vertical nodes. In any case, the g and r modes that 
are excited at an LR undergo wave channeling that is at 
least as strong as that of the f mode. 

In the case of young Jupiter in the solar nebula, GT80 
estimated that the wave-induced velocities are already 
mildly nonlinear (having Mach number about 0.3) at the 
LR. We therefore expect that shock damping would be 
important of a scale of order H from resonance, where 
wave channeling would occur. The extent of turbulent 
wave damping on that scale can be estimated by taking 
k ~ H~ x and a ~ 0.01. We then obtain the level of tur- 
bulent damping k\H <~ 0.1, and so we therefore expect 
that wave channeling strongly limits giant planet driven 
wave propagation in thermally stratified disks. 

8. SUMMARY AND DISCUSSION 
8.1. Summary of results 

We have analyzed the 3D response to tidal forcing of 
a gaseous disk having a vertical temperature variation. 
The unperturbed vertical disk structure was modeled as 
locally polytropic. We have considered effects at Lindblad 
resonances (LRs), subject to the thin disk approximation 
H <C r/m, for which the LRs provide horizontal forcing 
that is independent of height. In the standard vertically 
isothermal case considered to date, the horizontally prop- 
agating 2D mode (sound wave) carries all the resonantly 
generated angular momentum, but only if the thermody- 
namic response of the disk is also isothermal. 

In the vertically polytropic case, the f mode has a role 
equivalent to that of the 2D mode, although r and g 
modes are also launched at LRs (see Figure 1). The 
mode launched at vertical resonances (Lubow 1981) can 
be shown to be the p° mode. We have found that nearly 
all the torque exerted at an LR is carried by the f mode 
(see Table 1). Near resonance, the f mode behaves in a 
similar manner to the 2D mode in the isothermal case, 
in that it occupies the full vertical extent of the disk and 
behaves somewhat compressibly. The f mode is almost 
two-dimensional near resonance in that the vertical veloc- 
ities are smaller than the horizontal by a factor of order 
(iJ/r) 1 / 3 . However, away from resonance, the behaviors 
of the f and 2D modes differ radically (see §6.2). At a 
distance ~ r^/m from the resonance, the f mode begins 
to concentrate its energy near the surface of the disk, a 
process which we have called wave channeling. In this 
regime, the f mode behaves like a surface gravity mode 
(see Plates 1 and 2). In most cases, the wave amplitude 
increases by many orders of magnitude with distance from 
resonance, so that it is very likely that shocks will develop 
which would damp the wave. An important exception to 
this last statement occurs for the to = 2 ILR. The wave 
generated there undergoes relatively mild wave channeling 
and consequently a relatively mild increase in nonlinearity 
(by about a factor of 5). The amount of increase depends 
on the density distribution in the disk and could even re- 
sult in a decrease in dimensionless amplitude near the disk 
center. Dissipation by turbulent viscosity is likely impor- 
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tant in some cases, such as CV disks. On the other hand, 
the wave generated at the m = 2 OLR does exhibit strong 
wave channeling. These results indicate that circumbinary 
disks, as found around young binaries, and protoplanetary 
disks perturbed by planets, are subject to strong effects of 
wave channeling (see §7). 

8.2. Discussion 

The results in this paper differ from the previously ac- 
cepted picture for wave propagation in thermally stratified 
disks (e.g. LPS). The standard expectation was that a 
wave launched at a Lindblad resonance would begin prop- 
agating horizontally, but the wavefront would be rapidly 
refracted upwards as a result of the decrease in sound 
speed with increasing height above the mid-plane. Af- 
ter advancing a distance comparable to the disk thickness 
H, the wavefront would be substantially tilted upwards. 
The wave would then propagate vertically into the atmo- 
sphere of the disk, where it would shock. This model is 
based on the idea that the launched wave is a pressure 
wave or p mode in a high-frequency acoustic limit, so that 
one can consider the wavefront to propagate at the local 
sound speed without being affected by inertial or buoyancy 
forces. 

Our results provide a different picture. We have found 
that the wave cannot be considered to propagate verti- 
cally, since it is in fact a vertically evanescent f mode. In 
this view, taken in LP and KP, the disk behaves like a 
waveguide in which the wave is vertically confined. How- 
ever, somewhat similar to the standard picture, the wave 
energy does rise to the surface of the disk as the wave 
propagates away from the resonance. But this wave chan- 
neling process is effective over a distance of order r^/m, 
where tl is the radius of the Lindblad resonance and m 
is the azimuthal wavenumber. It does not depend on the 
disk thickness, provided that H/r <C 1. 

This can be understood roughly as follows. The wave 
channeling occurs because of the radial variation of the di- 
mensionless intrinsic frequency of the wave, (lo — mQ)/tt. 
This is equal to ±k/0 at the Lindblad resonance, but 
increases rapidly in magnitude over a distance of order 
tl/to. The wave therefore proceeds rapidly along the 
f-mode branch of the dispersion relation. In the high- 
frequency limit, as described in §6.1, the mode behaves 
like a surface gravity wave and is confined near the sur- 
face of the disk. (This would be true even in an incom- 
pressible disk in which refraction cannot operate.) The 
surface-gravity-mode dispersion relation is approximately 
(u> — mil) 2 g s k, where the surface gravity g s = fl 2 H and 
k is the radial wavenumber. It can be seen (see §6.1) that 
the mode is confined to the disk surface in a layer of verti- 
cal thickness 6 ~ k~ l . It then follows from the dispersion 
relation that (u> — mfl) 2 /f2 2 ~ H/5. As a result, the wave 
becomes confined (channeled) near the disk surface (i.e., 
8/ H < 1/2) over a radial distance from resonance of order 
r-^/m. 

Another issue concerns the boundary conditions. In 
LPS, the anticipated shocks in the disk atmosphere were 
represented in their simulation by an upper boundary con- 
dition that included some amount of dissipation. In our 
analysis of a polytropic disk, we have applied a bound- 



ary condition at the disk surface that acts to reflect waves 
rather than absorb energy. Even in a purely isothermal 
disk without a definite surface, as demonstrated by LP, 
all the modes are confined vertically and do not propagate 
vertically to infinity. This is a consequence of the increase 
in vertical gravity with height. 

Our treatment of the boundary is valid in the limit of 
high optical depth, where the atmosphere occupies a neg- 
ligible mass. In particular, the effects of atmosphere are 
not strongly felt by the wave until fc-ff atmos > 1, with 
ff a t m os being the density scale height at the base of the 
atmosphere. For a highly optically thick disk, H is much 
greater than £f a tmos, and so kH » 1. When this condition 
occurs, the wave has been strongly channeled before much 
wave energy enters the atmosphere. We plan to extend our 
current analysis to include the effects of a disk atmosphere 
in a future paper. 

The asymmetry in wave properties between the inner 
and outer LRs is most pronounced for low- to cases. For 
fixed to, the level of nonlinearity for ILR waves is lower 
than for OLR waves (see §6.2.2). One reason is that the 
density increases as the ILR wave propagates inward, while 
the density decreases as the OLR wave propagates out- 
wards (see eq. f57j). The other cause of the asymmetry 
is that the f-mode dispersion relation is followed to higher 
wavenumbers in the exterior of the OLR than in the in- 
terior of the ILR, for low m (see Figure 3). Away from 
the LR, the f mode is vertically confined to a surface layer 
of thickness ~ fc _1 . Due to wave action conservation, the 
confinement acts to increase the wave amplitude, causing 
substantial asymmetry in the case m = 2. In the high- 
m case, the behavior of the wavenumber becomes similar 
near the inner and outer LRs, and consequently the non- 
linearity becomes important at comparable distances from 
the respective resonances. 

A variant on the vertically polytropic model occurs for 
layered disks (Gammie 1996). In this model, some regions 
of protostellar or protoplanetary disks have turbulence re- 
stricted to the upper layers of the disk. The resulting 
vertical temperature structure is isothermal in the non- 
turbulent layer near mid-plane and the temperature de- 
clines with height in the outer turbulent layers. A wave 
launched from an LR in such a disk would probably re- 
quire a longer distance from resonance to undergo wave 
channeling. 

Given the large increases in wave velocities associated 
with wave channeling (see Plates 1 and 2), we expect that 
nonlinearities will cause the waves to shock and deposit 
angular momentum preferentially near the surface of the 
disk. We have described where nonlinearities are expected 
to become important in §6.2, based on a quasilinear esti- 
mate. Little work has been carried out to investigate the 
nonlinear behavior of f modes or surface gravity waves in a 
compressible fluid. It would be useful to conduct nonlinear 
simulations of the waves in such a disk. 
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Fig. 1. — Propagation of non-axisymmctric waves in a Kcplcrian disk. The five radii marked are the corotation radius ('CR'), the inner and 
outer Lindblad radii ('ILR', 'OLR'), and the inner and outer (vertical) resonances of the first p mode of even symmetry ('IVR', 'OVR'). The 
diagram indicates where the different classes of modes are excited and where they propagate. (The parameters used here are m = 2, s = 3, 
and 7 = 5/3.) 
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Fig. 3. — Variation of the dimensionless radial wavenumber with radius for the f° mode in a Keplerian, polytropic disk with s = 3 and 
7 = 5/3. The nine curves plotted correspond to m = 2, 3, . . . , 10, from bottom to top. 




Plate 1. (See Section 6.3 for discussion.) 
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Plate 2. (See Section 6.3 for discussion.) 



